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Abstract. 

This contribution focuses on the results of statistical model calculations 
at RHIC energies, including recently available experimental data. Previous 
calculations of particle yield ratios showed good agreement with measurements 
at SPS and lower energies, suggesting that the composite system possesses a high 
degree of chemical equilibrium at freeze-out. The effect of feeddown contamination 
on the model parameters is discussed, and the sensitivity of individual ratios to 
the model parameters (T, fig) is illustrated. 


1. Introduction 

Over the last few years, the study of thermalization in heavy-ion collisions using 
statistical models has been extended to SPS [lj, |], |jj and RHIC [| energies, allowing 
the so-called “freeze-out curve” || of hadronic matter to be extended close to the 
fig = 0 axis of the QCD phase diagram. This is an important development in the 
search for the quark gluon plasma because it allows a more direct comparison between 
lattice QCD calculations of the phase boundary temperature Jt|, performed mostly 
at zero baryon density, and experiment-based methods for measuring the equilibriun 
temperature on the hadronic side of the phase boundary (if the system equilibrates). 

The updated phase diagram presented at this conference by Braun-Munzinger 
Q suggests a scenario in which the final-state composition of the system, determined 
at chemical freeze-out, relies less on hadronic rescattering processes || as the beam 
energy is increased. This effect arises due to two reasons: stronger transverse flow 
causes faster expansion of the system, allowing less time for inelastic interactions; and 
the hadronic phase is bounded above by the temperature of the phase boundary itself, 
meaning that the system cools very little between hadronization and freeze-out. 

In our previous paper [|J , we showed that particle ratios at midrapidity measured 
during the RHIC run at yfs = 130 GeV are consistent with a system in chemical 
equilibrium, as described by a statistical model. The best-fit model parameters 
(T = 174 ± 7 MeV, \xb = 46 ± 6 MeV) differ by a factor 5 in baryon chemical potential 
/zb from the fit at maximum SPS energy Q, while the temperature changes very little 
(SPS: T = 168±3 MeV, /zb = 266±5 MeV). The small change in temperature, together 
with the good agreement with lattice QCD predictions of T c |u| (T c = 170—190 MeV), 
suggests that the system is nearly borne into chemical equilibrium out of the phase 
transition 0 0- We also presented predictions for particle ratios at full RHIC 
energy y/s = 200 GeV. 

In this contribution, the model comparison for the RHIC run at y/s = 130 GeV is 
extended to include additional ratio measurements which have been recently published 
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or presented, including at this conference. The treatment and effect of feeddown 
correction to particle densities in the model is discussed. Finally, the question of 
which ratios drive the model fit to the data is investigated. It is worth noting that the 
model parameter values (T, /jlb) obtained here for the best fit to the data have not 
changed from the fit performed with the data available at the time of our publication. 

2. Overview of experimental results 

The primary tool for investigating chemical equilibration in heavy-ion collisions is the 
experimental measurement of particle production, and in particular ratios of particle 
species. All four RHIC experiments have published or presented preliminary results 
on particle ratios. Some experimental observations are [fliif : 

• +/- ratios are approaching unity at midrapidity, indicating that quark coalescence 
and pair creation processes play a greater role in final-state particle yields. 
However, the midrapidity region is not yet baryon free at yfs = 130. 

• No large trends have been reported in ratios as a function of pt in the relevant p t 
range. 

• p/p increases from « 0.4 at y = 2 to « 0.65 at midrapidity, with a plateau width 
of 2 rapidity units. tt~/ n + shows no rapidity dependence over the same region. 

Since the submission of our paper Q, several additional measurements of ratios 
have been presented for the RHIC 2000 data. The (f/h~ ratio was measured by 
the STAR Collaboratio n fb§ to be 0.021 ± 0.001 sta t ± 0.004 sys . The PHENIX 
Collaboration presented |15[ preliminary results for K~ /7r — at this conference (most 
central: K~/n~ = 0.15 ± 0.01 sta t ± 0.03 sys ). 

Also at this conference, an updated preliminary value for K*°/h~ was presented 
Jill by the STAR Collaboration (|A'*° /h~ = 0.032 ± 0.003 sta t ± 0.008 sys ). The new 
value differs from the previous value Jl7| by nearly a factor 2 due to improved statistics 
and a different treatment of the residual background. 

3. Statistical model calculations 

The statistical model employed in the present calculations was first described in Q. 
In the model, particle densities are calculated as a grand canonical ensemble with the 
requirement that baryon, strangeness and isospin quantum numbers are conserved. 
All known particles and resonances up to 2 GeV are included in the calculation, and 
particles decay according to their branching ratios. 

The best fit model parameters (T, ps) are determined by comparing model ratios 
to experimental values and finding the minimum y. A y 2 /d.o.f. ~ 1 indicates a good 
fit. When both statistical and systematic errors have been presented for experimental 
values, the errors are added in quadrature. 

In this work, the best fit parameters are determined to be T = 174 ± 7 MeV, 
Pb = 46±6 MeV, with a y 2 /d.o.f. ~ 0.35. The parameters are the same as those in our 
paper [Q, even though a new ratio is included and another significantly revised. The 
low y; 2 /d.o.f. can be attributed to the large systematic errors in the first measurements 
at RHIC; if the fit is performed with only statistical errors in the experimental values, 
X 2 /d.o.f. increases to ~ 1.6. The true x 2 /d.o.f. lies between these two values. 
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Figure 1. Comparison between experimental data in Au+Au collisions at 
y/s = 130 GeV and statistical model calculations of particle ratios. The 
parameters of the model are T = 174 Mev, fiB = 46 MeV, and were determined 
by fi ndinp- the n ara. meter s et with mi nimum y 2 . Experimental data are taken 
from 0, 0.00 0000 00 


Figure 1 shows the model values for T = 174 MeV, /is = 46 MeV, and the 
available data for the four experiments at RHIC. Note that many of the ratios are 
preliminary. Model predictions for ratios involving h~ in the denominator have 
increased slightly due to the corrected treatment of negative hadrons. Now, as in the 
experimental values, only the primary negative hadrons (p, ir~ and K~) are included. 
In addition to the new and updated ratios mentioned in section 0 experimental values 
are taken from [0[ [0, 0, 21, 0 |0; see the table of compiled ratios in our paper 

for specific references. For the model calculations, a feeddown correction factor of 
0.5 was applied for weak decays; see section 3T for a discussion of this factor. 

Two observations can be made regarding Figure 1. First, the experimental value 
for \K*°\/h~ now lies under the model value but within error bars. The previous 
underprediction of \K*°\/h~ was a hot topic of discussion during this conference. 
Second, the good agreement between model and experiment for 4>/h~ is a success of 
the statistical model; the model value can be treated as a prediction, since the model 
parameters here have not changed from our paper Q| . 

It should also be noted that preliminary values were shown |0 at this conference 
by the STAR collaboration for A/n~ and S - /7r“. The A/ir~ measurement agrees 
well with the prediction from present model calculations (A/-7T - : expt = 0.066 ± 
0.004 sta t ± 0.005 sys ; model = 0.059). However, S _ /7r _ is underpredicted by a factor 
of 2 (S - /7r“ : expt = 0.014 ± 0.001 sta t ± 0.003 sys ; model = 0.072). If the experimental 
value for S - / n~ holds, the large disagreement with the statistical model would be an 
exciting result. However, we stress that the experimental results were presented as 
preliminary. 
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Figure 2. Minimum \ 2 points for three values of the feeddown contamination 
factor. Surrounding each minimum are lines of constant . A contamination 
factor of 1 indicates that no feeddown contributions have been removed from the 
“fed” particle multiplicities. 


3.1. Influence of feeddown correction 

As discussed in jjj primary particle multiplicities are affected by contributions from 
weak decays. In order to compare model calculations of particle production with 
experimental data, the “contamination” from these weak decays needs to be accounted 
for. This is done in the model by multiplying the branching ratios for weakly decaying 
particles (A, 2, fl, etc.) by a correction factor. 

Correction factors can be applied to particles individually in the model. However, 
in lieu of experimental values for these quantities, the default value used here and in 
our paper [Q is 0.5, which means e.g. that half of the protons arising from A decays 
are not removed from the proton multiplicity. 

Figure 2 shows the location of minimum yfl on the T — ps diagram for three 
values of the weak decay correction factor. The lines surrounding each minimum are 
contours of constant % 2 . The temperature decreases as the contamination from weak 
decays increases, while /j,b changes very little. The decrease in T can be explained 
by the direct relation between temperature and particle production in the statistical 
model. As the multiplicity of primary particles ( p , p , n, etc.) is increasingly 
contaminated, fewer true protons and pions need to be produced in order to agree 
with experimental values. 

Of course, it is not realistic that all feeddown correction factors have the same 
value. The result of Figure 2 merely demonstrates that weak decay contamination has 
a nontrivial effect on the extracted model parameters, particularly the temperature. 





Evidence for chemical equilibration at RHIC 
3.2. Ratios which constrain model parameters 


5 


In this section we investigate the contribution of individual ratios to the good 
agreement between the statistical model and experimental data. The question arises 
as to which ratios are strongly dependent on the model parameters (T, /is). A ratio’s 
sensitivity to these parameters comes from the mass and quantum numbers of the two 
particle species. This can be expressed in the Boltzmann approximation as: 

ni gi (mi\ 3 / 2 D c r 

~ ( I (uo-moMT' P’Brdi PS^i Pi Hi (1) 

ri 2 <72 Vm 2 / e ^ 2 rn2 >/ 1 

where gt is the spin-isospin degeneracy factor (< 7 , = (2 Si + 1)(2 R + 1 )), and B , S , and 
I represent the baryon, strangeness and isospin quantum numbers, respectively. The 
chemical potential p r is the sum of the individual chemical potentials in the system. 

Figure 3 shows the model calculation for each ratio included in Figure 1 as a 
function of ps for two values of the temperature. The range of values is chosen to 
include the fitted parameter values at y/s = 130 GeV. The p/p, A/A, and K~ /K + 
ratios are strongly dependent on fig. In the case of p/p, equation [I] becomes 

( 2 ) 


-2 ub/T 


The pb dependence of anti-baryon/baryon ratios is strong whenever p is much less 
than T ; therefore, the dependence weakens for strange and multistrange baryons due 
to the additional ps term in the chemical potential. 

Conversely, figure 4 shows the temperature dependence of each ratio for two 
values of ps- The p/ tt~ ratio exhibits the strongest dependence on the temperature. 
The strong T dependence of p/ tt~ is due to the large mass difference between the 
two species, which dominates for small values of p. This is also the primary reason 
why multi-strange baryon/meson ratios are sensitive to the temperature, which we 
demonstrated in our paper j|]. 

The S+/S - ratio’s dependence on temperature results from the large contribution 
of the ps term in the chemical potential (S 3 - = —2). Equation |l]is expressed as: 

'—2p B - 4 ps 




exp 


T 


(3) 


Note that strangeness conservation in the collision does not require that p$ is zero. 


4. Conclusions 

In this contribution, the work of our previous paper Q was extended to include 
new and updated ratios from the RHIC experiments. The best-fit parameters of 
the statistical model calculations do not change when including the new data. The 
strong agreement between model and data for the available ratios supports the notion 
of a system in chemical equilibrium at freeze-out. 

Also, the uncorrected feeddown contamination from weak decays has a small 
but relevant effect on the values of the extracted parameters. The model temperature 
decreases as feeddown contamination of primary particles (p, 7 r, etc.) increases because 
less primaries need to be produced by the model to account for the experimental value. 
Finally, ps was shown to be sensitive to anti-baryon/baryon ratios, primarily p/p, and 
the p/ tt~ ratio shows a strong temperature dependence. 

This work was done in collaboration with P. Braun-Munzinger, K. Redlich, and 
J. Stachel. The author thanks H. Van Hees and J. Knoll for fruitful discussions 
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regarding the statistical equations, and P. Braun-Munzinger for his careful reading of 
the manuscript. 
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Figure 3. Calculations with the statistical model showing the dependence on fiB 
for each ratio in Figure 1 for two different values of the temperature. Experimental 
values, compiled when multiple measurements are available, are represented by 
horizontal lines; shaded areas indicate the experimental error. 
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Figure 4. Calculations with the statistical model showing the temperature 
dependence of each ratio in Figure 1 for two different values of fig. Experimental 
values, compiled when multiple measurements are available, are represented by 
horizontal lines; shaded areas indicate the experimental error. 















































































